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Abstract. We consider the resolution of the Af = 2 supersymmetric KdV equation with 
a — —2 (SKdV a =—2) from the Hirota formalism. For the first time, a bilinear form of the 
SKdV a =-2 equation is constructed. We construct multisoliton solutions and rational similarity 
solutions. 



1. Introduction 

Let us recall that N = 2 supersymmetric (SUSY) extensions of the KdV equation have been 
largely studied in the past [TJ El [31 01 [5] in terms of integrability conditions and different types 
of solutions. Such extensions are given as a one-parameter (a € M) family of Grassmann-valued 
partial differential equations with one dependent variable A(x, t; 9i, 0%) where the independent 
variables are given as a set of even (commuting) space x and time t variables and a set of odd 
(anticommuting) variables 0\, 6 2 . The dependent variable A is supposed to be an even superfield 
and satisfies 

A t = (-A xx + (a + 2)AD 1 D 2 A + (a - l)(D r A){D 2 A) + aA 3 ) x , (1) 

where D\, D 2 are the covariant superderivatives given by 

D i = 9 i d x + de i , i = 1,2, (2) 

and such that D\ = D\ = d x . Note that since Df = d x (i = 1,2), Di can be viewed as the 
square root of d x . In |1J, it was shown that equation (TTJ) is integrable for the special values 
a = -2,1,4. 

Since the odd variables satisfy Q\ = Q\ = 9\ 9 2 + 6*2^1 = 0, the dependent variable A admits 
the following Taylor expansion 

A(x, t; 9 1 ,6 2 ) = u(x, t) + 0i£i(s, t) + 9 2 &(x, t) + O&vfa t), (3) 

where u and v are commuting (bosonic) complex valued functions and £1 and £2 ar e 
anticommuting (fermionic) complex valued functions. 

Solutions of equation (fTJ have been obtained using different methods. Among them, an 
original approach [2] has adapted the classical symmetry reduction method to the SUSY context. 



Indeed in this approach, starting with the SUSY equation ([I]), it has been reduced to a set of only 
one bosonic and one fermionic equations using invariance properties associated with invariant 
superalgebras. 

In the following, we will take a = — 2 in equation ([1]) where it is well-known that we get 
solitons solutions as travelling wave solutions [21 El H] and also rational similarity solutions as 
in the classical case [61 [7] . Such solutions will be generalized in the SUSY case using the Hirota 
formalism. Note that the Hirota formalism was adapted to the case a = 1 and a = 4 [5] but was 
still an open problem for the case a = —2. 

Let us write the set of equations we are working with in this case. The equation dTJ becomes: 



A t = {—A xx - 3(DiA)(D 2 A) - 2A%, (4) 
and admits the decomposition 

ut + (u xx + 2u 3 + 3^ 2 ) x = 0, (5) 

v t + (v xx + 3v 2 + 6u 2 v + 3u 2 x ) x - (3e 2 &ix + 12u&&)* = 0, (6) 

£i;t + + 3(u + 2u 2 )£i + 3u x &) x = 0, (7) 

+ + 3(« + 2u 2 )i 2 - 3u^i)x = 0. (8) 



It is interesting to note that in the one-fermionic parameter bosonization [S] which consists 
of writing the fermionic fields in Q as £j = C/i = 1>2), with /j a bosonic complex function 
and £ an odd parameter such that Q 2 = 0, equations (0) and © decouples into purely bosonic 
equations. In particular, we see that (JSJ reduces to the classical mKdV equation. 

2. Hirota formalism and solutions for the case a = —2. 

The Hirota formalism is a well known process in the classical [TU1 [TT] and in the M = 1 SUSY 
cases [12 ^1131 [T4"l ll5]. This formalism has been used, in particular, to obtain iV-soliton solutions. 
We examine the possible generalizations in the M = 2 SUSY case following the idea given in [8]. 
We write the superfield ([3]) in the following form 

A(x,t;9 1 ,9 2 ) = u b (x,t;6 1 )+6 2 ^(M;^), (9) 

where u b ', are even and odd superfields, respectively. Comparing with ([3|) , we have 
u b = u + and ^* = £2 — ^iu.We thus introduce the following change of variable A = d x B, 
where 

B(x, t; e l7 e 2 ) = U b (x, t; 6{) + 6 2 n f (x, t; (10) 

to equalize the order of the equation with the number of appearence of the x derivative in the 
nonlinear terms. Equation (JH) thus becomes, after integrating once, 

B t = -B xxx - ZiDxB^DzBx) - 2B 3 X , (11) 

where the constant of integration is set to zero. Inserting the explicit form ()10p of B in equation 
(fTTj) , we get a set of coupled Af = 1 SUSY equations on U b and r]f : 

U b + U b xxx - HDiU h x + 2(U b ) 3 = 0, (12) 
4 + vLx ~ HDi4 ~ ?U h xx D<U b x + 6r,f(U b ) 2 = 0. (13) 



Now the standard strategy is to introduce a change of the dependent variables, U b and 77^, 
in such a way that equations (|12p and (|13p become quadratic. 



Since equation (fl~2|) is a modified mKdV equation, we use the change of variables: 



U b = -i\og( T -), »7 / = J Dilog(-) ) (14) 

where t\ = T\(x,t;6i) and r 2 = r 2 (x,i;#i) are both bosonic superfields. 
From dHD, we deduce that 

^ = iD lU b , (15) 

which imposes the constraint v = —iu x and £2 = on the original superfield A =([3j). The 
constraints on the field £2 and v are standard when dealing with chiral superfields [16J. A chiral 
superfield A satisfies the constraint D + A = where D + = \(D\ + iD 2 ). Here, we see that these 
constraints are imposed in order to produce a bilinear form. 

From (|14|). we see that ry = iD\U b which reduces the pair of equations (|12p and (|13|) to the 
mKdV equation 

U b + U b xxx + 2(U b f = 0. (16) 

Let us mention that (fl~6|) is still a SUSY equation and using the decomposition U b = U^(x,t) + 
6\(f{x,t), we get the classical mKdV equation for Uq and ip satisfies 

<Pt + Pxxx + Q(Uq. x ) 2 cp x = 0, (17) 

for which a particular solution is (p = £ Uq- x where C is an odd constant. Such a result is common 
when we deal with M = 1 SUSY KdV and mKdV equations [El Q3J H US] . 
We thus have a direct bilinearization of equations (fT2|) and (fl~3|) 



(^ + P, 3 )(r r r 2 )=0, (18) 

5P !r (TrT 2 ) = 0, (19) 

where 

SV^(t v t 2 ) = (D 01 - D @2 )(d Xl - d X2 ) n T l {x l ;Q 1 )T 2 {x2;Q 2 )\x=x 1 =x 2 ,e 1 =0 1 =e 2 , (20) 

is the super Hirota derivative and Do, i = d 0i + ®id Xi , i = 1,2. Equations (|18j) and (|19|) are a 
natural generalization to the SUSY case of the classical bilinear form of the mKdV equation. 

2.1. Super soliton solutions. 

The Hirota formalism helps us to recover the travelling wave solution or one super soliton solution 
but we also get the N super soliton solutions [HE]. Indeed, for the travelling wave solution, we 
take 

ri = l + aie* r 2 = l + &ie*, (21) 

where \& = (kx + u)t) + #i£ and a± and b\ are nonzero even parameters. Introducing t\ and r 2 
in equation (|19j) yields the following relation 

61 = -ax. (22) 

The dispersion relation 

u) + k 3 = 0, (23) 

is obtained from the equation (|18p from which we get ^ = k(x — K?t) + ^i^. It is expected since 
we are interested in travelling wave solutions. Since 

dF 

U b = F(t>) = F(y ) + e 1 (—\„ = * , (24) 



where $o = k(x — K 2 t), and u + = d x U b , we see that the fermionic solution £1 is essentially 
the derivative with respect to x of u. 

We easily recover a one super soliton choosing a\ = i and k = 1. In such a case, the it 
component of the original superfield A is the classical soliton solution of the mKdV equation 
given by u(x, t) = sech(x — t). 

Now let us exhibit new solutions known as N super soliton solutions [H [T21 13 HH [15] . 

Indeed, in order to find the 2 super soliton solution, we first take 

n = 1 + aie* 1 + a 2 e* 2 + aia 2 A 12 e* 1+ * 2 , (25) 
t 2 = 1 + ftie* 1 + bae* 2 + b^B^e* 1 ^ 2 , (26) 

where now = KjX + ^jt + ^iCi and and b{ are nonzero even parameters (i = 1,2). This type 
of solution as for effect of breacking translation invariance. Introducing these expressions in the 
mKdV bilinear form (|18p and (|19p . the first equation yields the expected dispersion relations: 



Wi + «| = 0, i = 1,2. (27) 
The second equation gives the following conditions: 

bi = -Oi, i = 1,2, (28) 

A u = fl u =(^*»y (29) 

\Kl + K 2 J 

and the standard condition [U [13] relating the anticommuting variables d and £2 given by 

(«l,«2 / 0) 

= «2Ci- ( 30 ) 

Finally, the r-functions are given by 

Tl = 1 + aie* 1 + aae* 2 + aia 2 A 12 e* 1+ * 2 , (31) 
r 2 = 1 - aie* 1 - a 2 e* 2 + a x a 2 A 12 e* 1+ * 2 . (32) 

We may enjoy the behavior of the u(x,t) part of the 2-soliton solution 

9(10cosh(|(t - Ax)) + 5cosh(i - x) + 8sinh(|(t - 4x)) + 4sinh(i - x)) 
~ 72 + 41cosh(|(3t - Ax)) + 81cosh(|(7t - Ax)) + 40sinh(|(3t - 4a;)) ' 

where the parameters k±, k 2 , oi, a 2 are chosen as K\ = 2k 2 = 1 and a\ = a 2 = i, so that 
^1 = (x — t) + #iCi and ^2 = — i) + #iC 2 . In Figure 1, we give the solution \u\ and the 
fermionic part — fi in £1 = £/i = (u x for t = —10,0, 10. 




Figure 1. The functions \u(x, t)\ (full curves) and —fi(x,t) (dashed curves) for t = 



-10,0,10. 



For the three super soliton solution, t\ and r 2 take the explicit form: 

n = 1 + axe* 1 + a 2 e* 2 + a 3 e^ + a x a 2 A 12 e^ + ^ + ai a 3 A 13 e^ + ^ 
+a 2 a 3 A 23 e^ + ^ + ai a 2 a 3 A 12 A 13 A 23 e^ + ^ + *\ 

T2 = 1 _ ^e* 1 - a 2 e* 2 - a 3 e* 3 + a^^e* 1 ^ 2 + aic^e* 1 ** 8 
+a 2 a 3 A 23 e^ + ^ - ai a 2 a 3 A 12 A 13 A 23 e^ + ^ + ^ , 



(34) 
(35) 



where 



i^t + fliCi (* = 1, 2, 3) and 01 Q3] 



A; 



^2 "I - fc'i 



(36) 
(37) 



for i,j = 1, 2, 3 (i j). We may enjoy the behavior of the u(x, t) part of the 3-soliton solution. 
In Figure 2, we show again \u\ and the fermionic part —fx in £i = £/i = (,u x for i = —15, 0, 15 . 
The parameters K\, k 2 , k 3 , a±, a 2 , a 3 are chosen as k± = W-n 2 = in 3 = 1 and a\ = a 2 = a 3 = i. 






Figure 2. The functions \u(x, t)\ (full curves) and —fi(x,t) (dashed curves) for t = —15,0, 15. 



We have thus shown how to construct the N super soliton solution (N = 1,2,3), using the 
bilinear form (|18|) and (|19p . by giving the explicit forms of the functions t\ and r 2 . The N 
super soliton solution (N > 3) is easily generalized using the constraints above [HHHHTTJ. We 
are presently working on a Mathematica program which generates the ./V soliton solution by 
constructing the r- functions T\ and r 2 for general N. 

2.2. Rational similarity solutions. 

The Hirota formalism can also be used to get rational similarity solutions [10] . We assume now a 
SUSY generalization where the dependent variables t\ and r 2 are polynomials in the independent 
variable z = t~z (x+9\(,) and time t. Note that z\q 1= q= t~zx which is the invariant independent 
variable of the classical mKdV equation under dilatation [7]. For example, we take 

Ti(x,t;6i)=n(z,t)=th, T 2 (x,t;6 1 ) = T 2 (z,t) = t(5 3 + 12), (38) 
which solves the bilinear system (|18|) and (|19p . We thus have the following form 

U b (x, t;6 1 )= U b {z, t) = i log (tl Z3+ ~ U ^j • (39) 

From a Taylor expansion around 6\ = we get, zq = xi~3, 

U b (z, t) = U b (z ,t) + H8 l( ;—U b (z, t)U =Zo (40) 



and since u b = u + = d x U b , we get u(x,t) = t 3d zo U b (zo,t) and £1 = £u x , as expected. 
The generalization to the infinite set of solutions given in [61 [T7] is direct and we get 

*(,,«> (41) 

where the functions Q n (z) are the Yablonskii-Vorob'ev polynomials [17\ [T8] define by the 
recurrence relation 

^Qn+iQn-i = zQ 2 n - 12 {Q n Q n ,- z - z - Q 2 n . s ) , (42) 

with Qo(^) = and Qi(-S) = 5. Note that the Yablonskii-Vorob'ev polynomials are used to 
construct the similarity solutions of the Painlev II equation [171 [TH] . 

The link with the Hirota formalism is obtained by letting t\ and T2 be polynomials in the 
independent variable z and the time variable t. In fact, we take the following series 

n(n+l) (n + l)(n + 2) 

n, n (z,t)=t 6 Q n (z), T 2 , n (z,t)=t 6 Q n+ i(z) (43) 



which lead to (|4T[) and solve the bilinear equations p8j) and (|T9|) . For example, we see that if 
n = 1, we recover the expression (|38p . 

From the Taylor expansion (|40p and (|4ip . we get the corresponding u n (x,t) solution of the 
bosonic superfield A, namely 

U n ( X ,t) =it-^\0g( Q ^ {Z " ) ) , C/,:„ <<),,!„. (44) 

dz V <3n(^o) / 

Figure 3 shows the behaviour of the imaginary part of u\ and f\-\ as a function of x and t, 

x 3 — 6t 

ui(x,t) = 2i—^ -. (45) 

v ; x(x 3 + 12t) v ' 




3. Conclusion. 

The generalization of the Hirota formalism to M = 2 supersymmetric equations has been used 
to solve the celebrated SUSY KdV equation for a = —2. We have been able, for the first time, 
to produce a bilinear form of the N = 2 SUSY KdV equation. In fact, the assumption made 
has produced a reduction of our equation to a N = 1 SUSY one for which the bilinearization 
is known. From the bilinear equations we have constructed multisoliton solutions and have 
produced for the first time a SUSY generalization of rational similarity solution. Indeed, we have 



generated an infinite set of rational solution using a SUSY version of the Yablonskii-Vorob'ev 
polynomials well known in the construction of similarity solution of the Painlev II equation. 

In future work, we want to use the symmetry reduction method adapted to the SUSY context 
to create links with the Hirota formalism. This paper suggest such links. Indeed the Yablonskii- 
Vorob'ev polynomials, which arise from the classical mKdV equation via dilatation invariance, 
were used in the Hirota formalism to generate an infinite set of solution of the TV = 2 SUSY 
KdV equation. 
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